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Technical Memorandum X-73349 


SPACELAB EXPERIMENT COMPUTER STUDY 
Volume rV: Spacelab User Cost Data ( Central Experiment Computer)' 


SECTION 1 Option lAl - Central .with Mini, Central Software 

Development by Central Group . 

1. 1 Costing Method 



ASSUMPTIONS 


o STILt facility hardware used to maximum extent in lieu of pur- 
chasing added hardware. 

o Mini’s required are procured such that they are built compatible 
with existing STIL interface hardware to accommodate real-time 
simulations for software checkout and verification. 

o DEP mini's are procured that have existing cross support soft- 
ware (assembler, link editor, Fortran Compiler) that will 
execute on the STIE Host (S360) 

o All DEP mini’s are. assumed to be procured from different 
vendors and are not software compatible 

o A centralized group of programmers will develop all software. 
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Cost Element 4. 1 Experiment Application Software Development 
Cost Factors 

1. Software Development 

2. Common Software 

3. Host Computer Time 

4. Simulation Computer Time 

5. Host Computer Time DEP Software 

- 6. Simulation Computer Time DEP Software 

7. Travel 

8. T raining 

(1) Software Development 

(A) = {{Number of Statements) {Cost/Statement)) Per 
Flight Per Yr. 

Number of Statements = GDC Estimate = {{New Flights) - 
{Previous Yr. Common)) Per Yr. . 

Cost/Statement = $45 {Requirements, Code, and Verification) 
$15 Per Statement is Applied to Integrated Verification 
{See 4. 3 {!)) 

{1) Total = {{#10) {$45)) Per Yr. 

{2) Common Software 


{A) {{Number of Statements) {Cost/Statement)) Per Yr. 
for 5 years Starting in FY81. 

Number of Statements = 10% of new Development per 
Flight Per Yr. Reducing 2% Per Yr. to 0% - Yearly 
Totals are Accumulative and Accumulated Total is 
Subtracted from Subsequent Years Required Development. 

Cost/Statement = $5 {Document and Place in Library) 
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(2) Total = ((#5) ($5)) Per Yr. 


(3) Host Computer Time 


(A) = ((Host Time) (Cost/Hr, )) Per New Flight Per 
Yr. 

Host Time = ((# Instructions) ((# Instructions /Module) 

(Hrs. /Module)) 

= (# of Modules) (Hrs. /ModTile) 

Hrs. /Module = 18 Compiles /Module @ 3 Mins. /Compile 
+ 9 Functional Simulations @ 12 Mins. / 
Run + 1 Data Reduction for 75% of 
Simulation Runs @ 10 Mins. Each 

= (54 Mins, ) + (108 Mins, ) + (70 Mins, ) 

= 3. 87 

Host Cost/Hr, = ((10% Maintenance) + (Operations) 

+ (Consumables)) -r (2080 Hrs. /Shift) 

= (2. 43) + (57, 69) + (63. 05) = $123. 22/Hr. 

# of Modules = #8 

(3) Total = (3.87) (#8) ($123. 22)/ Yr. 

(4) Simulation Computer Time 


(A) = (# of Hrs. /Module) (Cost/Hr. ) (# of Modules) 

# of Hrs. = 4 Simulations at 60 Min. /Simulation 
(Includes Set-Up, Runs, and Run Evaluations) 

Cost/ Hr. s: (10% (Maintenance (CDMS, CID, 
SIMULATION COMPUTER)) + (Consumables) 

+ (Operations)) -r 2080 Hrs, 
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■ = ($8) + (0) + {(2) ($40K))4 2080 

= ($8) + ($38.46) 

= $46. 46/Hr. 

(4) # of Modules = #8 

(4) Total = (($46,46) (#8) (4 Hrs, )) Per Yr. 

(5) Host Computer Time PEP Software 

(A) = ((Host Time/Module) (#Modules) (Cost/Hr. )) Per 
Flight/ Yr, 

Host Time = (# Modules) (Time/Module) 

Time/Module = 18 Compiles /Modtile @ 3 Mins. /Compile 
+ 3 Data Reduction Ruus/Module @ 10 Mins, /Run 

= 84 Mins. /Module 

Cost/Hr. = $123. 22 (Same as 4.1 (3)) 

if of Modules = #9 

(5) Total = ({1. 4 Hrs. ) (f/9) ($123. 22)) /Yr. 

(6) Simulation Computer Time DEP Software 

(A) =((# of Hrs. ) (Cost/Hr. ) (# of ModTiles)) Per Flight/Yr. 

# of Hrs. = 7 Simulations at 60 Mins. /Simulation 
(Includes Set-up, etc. ) 

* 

9 Functional Simulations = 3 Realtime Simxilations . 
Assume no Functional Simulators for DEP, 

Cost/Hr. = $46.46 (Same as 4.1 (4)) 

(6) Total = {(7 Hrs.) ($46.46) (#9)) Per Yr. 
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(7) Travel 


(A) = ((Number of Man Yrs. ) (Travel Cost/Man Yr. ))/ Yr. 

Nimiber of Man Yrs. = Number of Man Yrs. Required 
as Determined by Number of New Instructions. 

Travel/Man Yr. = Number of Trips ((Cost of Ticket) + 
(Number of Days) (Cost Day)) 

Number of Trips = 4/Man Yr. 

Cost of Ticket - $150 

Cost/Day = $12. 50 Per Diem + $12. 50 for Car = $25 

Number of Days = 365 

Travel = 4(($150) + (365) ($25)) 

= $600 + $9,125 = $9, 725/Man/ Yr. 

(7) Total = 0 for This Option 
(8) Training 

(A) = (Number of Programmers) (Cost/Programmer) 

Niimber of Programmers = Number of Programmers 
Required by Number of New Instructions. 

Cost of Programmer = Engineering- Estimate $500/Man 

(8) Total = 0 for This Option 
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Cost Element 4. 2 Experiment Application Software Maintenance 


Cost Factor_s 

1. Experiment Unique Software 

2. Experiment Common Software 

3. Host Computer Time 

4. Simulation Computer Time 

5. Host Computer Time DEP Software 

6. Simulation Computer Time DEP Software 

7. Travel 

(1) Experiment Unique Software 

(A) = {# of Statements) (Hate of Change) (Cost/Statnment) 
(Number of Statements) (Rate of Change) = //1 3 

Rate of Change = Engineering Estimate 13ased on Past 
Programs - $0% for 1st Re-Fly, 30% for 2nd Re- Fly, 

20% for 3rd Re-Fly, 10% for all Subsequent Re-Flys. 

Cost/Statement = $45 (See 4.1 for Rationale) 

(1) Total = ((#13) ($45)) Per Yr. 

(2) Experiment Common Software 


(A) = ((// of Statements) (Change Rato) (Cost/Statemont)) 
Per Yr. 

Rate of Change = Engineering Estimate = 1% 

Cost/Statement = $60 (Due to Verification Necessary for 
Multi-Use) 

(2) Total = ((//6) (1%) ($60)) Per Yr. 
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(3) Host Computer Time 

(A) = (#Modtiles) (#Hrs. /Module) (Cost/Hr. ) 

#Modviles = Line {#14) 

§ Hrs, /Module = 3. 87 (See 4,1 for Rationale) 

Cost/Hr. = $123,22 (See 4.1 for Rationale) 

(3) Total = (( #14) (3. 87) ($123.22)) Per Yr. 

(4) Simulation Computer Time 

(A) = {(#Modules) ( #of Hrs. /Module) (Cost/Hr, )) Per Yr. 

# Modules = Line (#14) 

# Hrs. /Module = 4 Hrs. (See 4,1 for Rationale) 

(4) Total = ((#14) (4) ($46.46)) Per Yr. 

(5) Host Computer Time DEP Software 

(A) ((Host Time/Module) (#Modules) (Cost/Hr. )) Per Yr. 
Host Time/Module = 1.4 Hrs. (See 4.1 for Rationale) 

# Modules = Line (#15) 

Cost/Hr. = $123. 22 (See 4.1 for Rationale) 

(5) Total = ((1.4) ( #15) ($123.22)) Per Yr. 

(6) Simulation Computer Time DEP Software 

(A) ((# of Hrs. ) (Cost/Hr.) (#Modules)) Per Yr. 

# of Hrs. = 7 (See 4.1 for Rationale) 

Cost/Hr. = $46,46 (See 4.1 for Rationale) 

(6) Total = ((7) ($46.46) (#15)) Per Yr. 
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(7) Travel 


(A) ((Number of Man Yrs. ) (Travel Cost/Man Yr. })/Yr. 

Number of Man Yrs. = Number of Man Yrs. Required as 
Determined by Number of Instructions to be Maintained. 

Travel = $9, 725 Man/Yr. 

(7) Total = 0 for This Option 
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Cost Element 4. 3 


EAS Software Integrated Verification 


Cost Factors 

1. Integrated Verification 

2 . Host Computer Time 

3. Simulation Computer Time 

4. Integrated Verification Simulation Software 

(1) Integrated Verification 

(A) = ((# of Modules) (Cost/Modtile)) Per Flight/ Yr. 

# of Modules = #21 

Cost/Module ^ {# of Statement) (Cost/Statement) 

# of Statements = 100/Module 

Cost/ Statement = $15 ($60 Estimated Total Cost Per 
Statement for Central Development Less $45 for 
Development) 

= $1, 5 00 /Module 

(1) Total = ((#21) ($1,500))/Yr. 

(2) Host Computer Time 

(A) = {# Hrs. ) (Cost/Hr, ) 

(2) Total = (0) For This Function (Has Been Included in 
Required Host Runs 4.1) 

(3) Simulation Hardware Time 

(A) = (# Hrs. /Module) (Cost/Hr. ) (# Module) 

# Hrs. /Module = 2 Simulations /Module at 60 Min. Per 
Simtalation (Includes Set-Up, etc.) 
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Cost/Hr. = $46.46 {See 4,1 (4) B) 

# Modules = #21 

(3) Total = ((2) (#21) ($46.46)) Per Yr, 

(4) Integrated Verification Simulation Software 
(A) = ('(# of Modules) (Cost/Module)) Per Yr. 

# of Modiiles = #29 

Cost/Module = $1, 500 (See 4. 3 (1) for Rationale 

(4) Total = ((#29) ($1,500)) Per Yr. 
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Cost Element 4. 4 


Prefltght Checkout Software Development 


Cost Factors 

1. Software Development 

2. Common Software 

3. Host Computer Time 

4. Simulation Computer Time 

5. Travel . 

(1) Software Development 

(A) = ((Number of HOL Statements) (Cost/Statement)) Yr, 

(1) Total = ((#22) ($30))/ Yr, 

(2) Common Software 

(2) Total = 0 for This Option 

(3) Host Computer Time 

Host Computer Time = ((Host Time)) (Cost/Hr, ) 

Kef. 4,1 (3) 

(A) = ((3,87) (# Modules) ($123, 22))/ Yr. 

(3) Total = {(3. 87) (#22 4- 100) ($123. 22)) Per Yr. 

(4) Simulation Computer Time 

(A) = (# Hrs. /Module) (Cost/Hr. ) (# Modules) 

Ref. 4.1 (4) 

(4) Total = (4) ($46. 46) (#22 100) Per Yr. 
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TABLE 3 


Buoys Used in Program LOCATE 


Buoy ID 


North Latitude- 


West Longitude 


41001 

35.0 

41002 

32.3 

41004 

32.6 

42001 

26.0 

42002 

26.0 

42003 

26.0- 

44001 

38.7 

44002 

40.1 

44003 

40.8 

44004 

39.0 

46001 

56.0 

46002 

42.5 

46003 

52.0 

46004 

51.0 

46005 

46.0 

46006 

41.0 

46007 

59.2 

46008 

57.1 

46009 

60.2 


72.0 
75.3 

78.7 

90.0 

93.5 

86.0 

73.6 

73.0 
68.5 

70.0 

148.0 

130.0 

156.0 

136.0 

131.0 

138.0 

152.7 
15T.7 

146.8 
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SMMR TB (K) 


Figure 13.10. Orbit 1212, 38“ N, Nominal and Interim 


SMMR 37.0 H TB VS LATITUDE 
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SMMR TB (K) 


Figure 13.9, Orbit 1212, 38® N, Nominal and Interim 


SmR 37.0 V TB VS LRT3TUDE 
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SMMR TB (K) 


Figure 13>8. Orbit 1212, 38® N, Nominal and Interim 

SmR 21.0 H TB VS LBT]TUDE 



lRTITUDE (DE6.) 
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(K) 


Figure 13.7. Orbit 1212, 38® N, Nominal and Interim 

SMMR 21.0 V TB VS LRT3TUDE 




SMMR TB (K) 


Figure 13.6. Orbit 1212, 38° N, Nominal and Interim 


SmR 18.0 H TB VS LRT3TUDE 
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SMMR TB <K) 


Figure 13.5* Orbit 1212, 38® N, Nominal and Interim 
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SMMR TB (K> 


Figure 13.4. Orbit 1212 , 38® N, Ncjminal and Interim 


SMMR 10.7 H TB VS LRTITUDE 



lRTITUDE (DEG.) 


129 



SMMR TB (K) 


Figure 13*3. Orbit 1212, 38° N, Nominal and Interim 

SmR 10.7 V TB vs LATITUDE 



LATITUDE (DEG.) 
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SMMR TB (K) 


Figure 13* 2* Orbit 1212, 38® N, Nominal and Interim' 


SWR 6,6 H TB VS LnilTUDE 
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SMMR 


Figure 13.1. Orbit 1212, 38® N, Nominal and Interim 


SMMR 6.6 V .TB VS LATITUDE 



34. 36. 3a. 40. 42. 


LATITUDE (DEG.) 
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SMMR TB (K) 


Figure 12,7. Orbit 1212, 38® N, Cross and Interim 


SriMR 21.0 V T5 VS LATITUDE 
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SMMR TB (K) 


Figtjre 12.6. Orbit 1212, 38° N, Cross and Interim 


ShMR ia,0 H TB VS LRTITUDE 
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SMNF^ ID no 


Figure 12.5. Orbit 1212, 
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r^l 


SMMR TB (K) 
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r c\j 


Figure 12.3. Orbit 1212, 38® H, Cross and Interim 


SmR 10.7 V T5 vs L(=^TITUDE 



latitude (DEG.) 
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SMMR TB (K) 


Figure 12.2. Orbit 1212, 38'" N, Cross aad Interim 


SKMR 6.6 H TB VS LfiTITUDE 
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SMMR TB (K) 


Figure 12.1. Orbit 1212, 38® N, Cross and Interim 


SMMR 6.6 V T5 VS LRTITUDE 
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SMMR TB (K) 


Figure Orbit 1212, 38® N, Box and Interim 

SMMR 37.0 H TB VS LATITUDE 



latitude (DEG.) 
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SMMR 


Figure 11.9. Orbit 1212, 38° H, Box and Interim 

SMMR 37.0 V TB VS LRTITUDE 



LRTITUDE (DEG.) 


114 




SMMR TB (K) 


Figure 11.8. Orbit 1212a 38® N, Box and Interim 


SMMR 21.0 H T5 VS LATITUDE 



latitude (DEG.) 
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(K) 


Tigure 11 * 7 * Orbit 1212, 38® N, Box and Interim 




SMMR TB (K) 


Figure 11. 6* Orbit. 1212, 38® N, Box and Interim 


SMMR 18.0 H TB VS LBTITUDE 
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SMMR 


Figure 11*5. Orbit 1212 > 38* N, Box and Interim 

SMMR 16.0 V TB VS LRTlTUdE 



56. 53. 4-0- 42. 

lPTITUDE (DEG.) 
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SMMR TB (K1 


11.4. Orbit 1212j 38* N, Box and Interiai 


SmR 10.7 H TB vs LRTITUDE 
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(K) 


Figure 11,3, Orbit 1212, 38® N, Box and Interim 

SMMR iO.7 V TB VS LRTITUDE 




SMMR TB (K) 


Figure 11.2. Orbit 1212, 38® N, Box and Interim 


SMMR 6.6 H TB VS LRTITUDE 



lRTITUDE (DEG.) 


107 



SMMR TB (K) 


Figure 11.1. Orbit 1212, 38° N, Box and Interim 

SMMR 6.6 V TB VS LRTITUDE 
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SMMR TB IK\ 


Figure 10.10 Orbit 1212, 55° N, Noniinal and Interim 


SWR 37.0 H TB VS LATITUDE 



LATITUDE (DEG.) 
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rbit 1212, 55® N, Nominal and Interim 

IR 37.0 V TB VS LATITUDE 






ShMR TB n<l 


Figure 10.8. Orbit 1212, 55° N, Nominal and Interim 


SMMR 21.0 H T5 VS LATITUDE 
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SMMR TB (Ki. 


Figure 10.7. Orbit 1212, 55“ N, Nominal and Int 
SMMR 21.0 V TB VS LRTITUOl 
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SMMR TB IK) 


Figure 10.6. Orbit 1212, 55® N, Nominal and Interim 


SMMR la.Q H T5 VS LATITUDE 



latitude (DEG.) 
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. in 








SMMR TB IK) 
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. to 



cbit 1212, 55* N, Nominal and Interim 

IR 6.6 V TB VS LATITUDE 




SMMR TB IK\ 


Figure 9.9. Orbit 1212, 55“ N, Cross and Interim 


SmR 37.0 V TS VS LR7ITUDE 
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SMMR TB (K\ 


Figure 9.8, Orbit 1212, 55® N, Cross and Interim 


SmR 21.Q H T5 VS LATITUDE 
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SnhR TB (Kl 


Figure 9.6. Orbit 1212, 55° N, Cross and Interim 


SMriR la.O H T6 vs LATITUDE 



91 






SMMR TB lk\ 


Figure 9.4. Orbit 1212, 55" N, Cross and Interim 


SHMR 10 7 H TS VS LATITUDE 
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(K1 


Figure 9,3. Orbit 1212, 55® N, Cross and Interim 




SMMR TB (kl 


Figure 9.2. Orbit. 1212, 55® N, Cross and Interim 


SMiiR 6.6 H T5 VS LRTITUDE 
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SMMR 
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SMMR 


Figure 8.9. Orbit 1212 j 55° N, Box and Interim 


. SMMR 37.0 V T5 VS LflTITUDE 
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Figure 8.7. Orbit 1212, 55® N, Box and Interim 


SWR 21 0, V T5 VS latitude 



SmR 1 6 IK) 


Figure 8.6. Orbit 1212, ! 


F=tiT*i 



L.ciTITUDE (DEG.) 
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SMMR ’lb IK) 


Figure 8*5. Orbit 1212, 1 


SMnR 18.0 V TB VS LATITUDE 



latitude (DEG.) 
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8i 




SMMR Tb no 


Figure 8.4. Orbit 1212, 55° N, Box and Interim 


SMMR 10 7 H TB VS LPTITUDS 



latitude (DEG.) 
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SMMR TB no 


Figure 8*3, Orbit 1212, 55® N, Box and Interim 


SM.MR 10.7 V T5 VS LRTITUDE 
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SMMR ’tri no 


Figure 8.2. Orbit 1212, 55® N, Box and Interim 


SmR 6.6 H TB VS.LR7ITUDE 
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ShMR (K) 


Figure 8,1. Orbit 1212, 55® K, Box and, Interim 


SMMR 6.6 V TB VS LRTITUDE 
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Figure 7.5. Orbit 331, Nominal Mode 


SMMR 37.0 GHZ T5 CROSS TRRCK GRADIENT VS LATITUDE 



0. 10. 20 . ' ' 30. 40. 

LATITUDE (DEG) 
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Figure 7.4. Orbit 331, Nominal Mode 


SMf^R 21 0 GHZ T5 CROSS TRRCK GRPDIENT VS LATITUDE 



LATITUDE (DEG) 
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Figure 7.3. Orbit 331, Nominal Mode 


SMMR 13.0 GHZ TB CROSS TRRCK GRRDIENT VS LATITUDE 



LATITUDE (DEG) 
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Figure 7.2. Orbit 331, Nominal Mode 

SnnR 10.69 GHZ TB CROSS TRRCK GRADIENT VS LfiTiTUDE 
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Figure 7.1. Orbit 331, Nominal Mode' 


SMMR 6.6 GHt T5 CROSS TRRCK GRRDIENT VS LRTITUDE 



LATITUDE (DEG) . 
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Pigure 6.5* Orbit, 331, Cross Mode 


SMMR 37.0 GH2 T5 'CROSS TRRCK GRADIENT VS LATITUDE 
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Figure 6.4. Orbit 331, Cross Mode 


SM^lR 21.0 GHZ T5 CROSS TRACK GRADIENT VS LATITUDE 



LATITUDE (DEG) 
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Figure 6.3. Orbit 331, Cross Mode 


SMM{R 18.0 GHZ T5 CROSS TRACK GRADIENT VS LATITUDE 
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Figure 6.2. Orbit 331, Cross Mode 


SMHR 10.69 GHZ T5 CROSS TRACK GRADIENT VS LATITUDE 



LffTITUDE (DEG) 
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Figure 6.1. Orbit 331, Cross Mode 


SMMR 6.6 GHZ T5 CROSS TRACK GRADIENT VS LATITUDE 
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Figure 5.5. Orbit 331, Box Mode 

SMMR 37.0 GHZ T5 CROSS TRACK GRADIENT VS LATITUDE 
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Figure 5*4* Orbit 331, Box Mode 


SliMR 21.0 GHZ T5 CROSS TRfiCK GRRDJENT VS LATITUDE 



LATITUDE (DEG) 


64 







Figure 5.3. Orbit 331, Box Mode 


SMMR 18.0 G.HZ T5 CROSS TRACK GRADIENT VS LATITUDE 
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Figtire 5.2. Orbit 331j Box Mode 
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TB CROSS-TRFICK GRRD1E:NT- (K/CEXL’) 


Pigure 5.1* Orbit 331, Box Mode 


SMMR 6.6 GHZ TB CROSS TRfiCK GRADIENT VS LATITUDE 
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Figure 4.5. Orbit 331, Interim Mode 
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Figure 4.4. Orbit 331, Interim Mode 


SMMR 21.0 GHZ TB CROSS TRRCK GRADIENT VS .LATITUDE 



LATITUDE (DEG) 


59 



Figure 4.3, Orbit 331, Interim Mode 


SMMR 18.0 GHZ T5 CROSS TRPCK GRRDIENT VS LATITUDE 



LATITUDE (DEG) 
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TB CROSS-TRRCK GRODIETMT (k/CELL) 


Figure 4.2. Orbit 331, Interim Mode 


SMMR 10.69 GHZ T5 CROSS TRfiCK GRADIENT VS LfiTITUDE 



LATITUDE (DEG) 
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Figure 4,1. Orbit 331, Interim Mode 
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SMMR TB iK) 


Figure. 3. 10. Orbit 331, Nominal and Interim 

SMMR 37.0 H T5 VS LATITUDE 
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l\<\ 


Figure 3.8. Orbit 331, Nominal and Interim 

SMMR 21.0 H T5 VS LPTITUDE 






SMMR TB (Kl 


Figure 3.6. Orbit 331, Nominal and Interim 

SMMR 16.0 H T5 VS LRTITUDE 
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Figur 
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SMMR TB n<) 


Figure Orbit 331, Nominal and Interim 


SMMR 6.6 V T5 VS LfiTlTUDE 
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SMMR TB IK) 


Figure 2.7. Orbit 331, Cross and Interim 


SMMR 21.0 V TS VS LATITUDE 
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SMMR hi IK) 


Figur 
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SMMR ‘I 


Figure 2.2. Orbit 331, Cross and Interim 


SmR 6.6 H T3 VS •LATITUDE 



37 



SHMR 


Figure 2.1. Orbit 331, Cross and Interim 


SMMR 6.6 V T5 VS LRTITUDE 



LATITUDE (DEG.) 
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SMMR 
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SMMR IB no 


Figure 1.9. Orbit 331 


, Box and Interim 

SMMR 37.0 V T5 VS LPTITUDE 
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SMMR TB M<^ 


Figure 1.8. Orbit 331, Box and Interim 


SmR 21.0 H T5 VS LfiTITUDE 
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SMMR TB no 


Figure 1.6. Orbit 331, Box and Interim 

SliHR ia.O H T5 vs LATITUDE 
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SMMR TB IK\ 


Figure 1.5. Orbit 331 


, Box and Interim 


SMMR IB.O V T5 VS LATITUDE 
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SMMR Tli (\<\ 


Figure 1.4. Orbit 331, Box and Interim 


SMMR 10.7 H T5 VS LRTITUDE 
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SMMR TB no 


Figure 1.2. Orbit 331, Box and Interim 


SmR 6.6 H TB VS LATITUDE 



LATITUDE (DEG.) 
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SMMR TB MO 


Figure 1*1. Orbit 331, Box and Interim 


SMMR 6.6 V T5 VS LATITUDE 



26 




Table 2* Land Sidelobe Effects Orbit 1212, 51*N to 59*^N, Colntm 1 


INTERIH BOX INTERIM - CROSS 'IStERIM - NOMI NAL 



Open 


Near 

Open 


Near 

Open 

Near 

Channel 

Ocean 

N* 

Landt 

Ocean 

N* 

LandT 

Ocean 

Landt 

6.6V 

4.50 

7 

6.50 

4.50 

7 

8.50 

4.50 

15 


6.6H • 

3.25 

7 

5.00 

3.25 

7 

6.75 

2.75- 

18 

** 

10.7V 

3.00 

8 

5.00 

3.00 

8 

7.00 

3.75 

18 

** 

10. 7H 

4.00 

8 

6.00 

4.00 

10 

9.00*** 

4.00 

>20 

** 

18V 

3.00 

4 

4.00 

3.50 

4 

5.50*** 

5. do' 

9 


18H 

4.50 

6 

5.50 

4.50 

6 

8.00 

5.00 

>28 


21V 

3.50 

0 

3.50 

3.50 

0 

3’. 50 

3.50 

14 


21H 

3.75 

5 

5.00 

3.75 

. 5 

5.50 

4.00 

22 

** 

37V 

3.00 

0 

3.00 

3.00 

3 

4.00 

4.00 

24 


37H 

4.00 

4 

5.00 

4.00 

4 

7.00 

4.00 

>20 

** 


t Near Land refers to cell adjacent to land. 

* N is the number of cells from land at which the difference first equals 
the *'open ocean difference". 

** The near land values for the nominal mode exhibit ringing effects, 

*** CROSS mode appears to be over correcting for nearby land. 

(Box and Cross appear to be identical except for the cell adjacent to land). 
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Table 1.3. Orbit 331, 33®N to 43^N Interim Minns Nominal Mode (Continued) 


channel column 


37i0 V 1 

2 
3 
H 

5 

6 

7 

8 

9 

10 
I ] 
12 
13 
H 

15 

16 
17 
IS 

19 

20 
21 
22 

37.0 H I 

2 
3 

5 

6 

7 

8 
9 

10 
i 1 
12 
13 

n 

15 

16 
17 
16 

19 

20 
21 
22 


mean 

STD. DEV 

3. so 

2.36 

3 . 6 9 

2.81 

3.61 

3.79 

2.73 

3.10 

2. *13 

2.27 

2*2^ 

1.53 

2.2q 

1.7m 

2.23 

2.18 

1 .79 

2. ‘19 

2.6M 

2.91 

2.13 

2.75 

1 .97 

1 .50 

2.M6 

2.7M 

2.03 

2.71 

1.96 

2.^7 

2.32 

2.11 

1 .71 

1.93 

2*11 

1 .*16 

2.02 

1 .85 

2.16 

2.27 

1.96 

1.39 

2.65 

3.2M 

3.3^ 

2-.63 

3t0i 

M.75 

2.85 


2.93 

M.08. 


1 i93 

3.07 

2.9S 

3. SO 

2.77 

2.55 

3.68 

3.20 

3.28 

3.06 

2. IM 

2.57 

3.55 

2.90 

2.76 

2.92 

3.55 

2.87 

3.05 

3.01 

2.0M 

3.08 

3.85 

2.58 

2.2m 

3.56 

1 .99 

3.H8 

5* 15 

3.01 

T 1 .90 

H.63 

6.33 

2t5H 

13.7m 
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Table 1 


Orbit 331, 

33°N to 43®!^ 

Interim Minus Nominal Mode (Continued) 

channel 

COLUMN 

MEAN 

STD. DEV« 

18«Q H 

1 

*» • *13 

1*16 

' 

2 

3.38 

1.02 


3 

3.09 

1.25 


M 

3.08 

1.29 


5 

2.98 

1 .Ms 


6 

2.86 

1.32 


7 

3.06 

1.55 

• 

8 

3.15 

1.22 


9 

3.11 

. 06 


10 

3.35 

.83 


U 

3.93 

1*12 

21 *0 V 

1 

^ .99 

1.62 


2 

3.60 

♦ 96 


3 

3.70 

.83 


H 

3.69 

.89 


S 

3.50 

1 .06 


6 

3;so 

.88 


7 

3. M3 

1 .mm 


S 

3.86 

.92 


•9 

3.55 

.63 


10 

3.66 

.77 

. 

11 

M.37 

.95 

21 .Q H 

1 

M,25 

1.57 


2 

3.3M 

1.11 


3 


1 .58 


9 

3.70 

1.37 


5 

3.3M 

I.IM 


6 

3.05 

1 .5-8 


7 

3.19 

1.19 


B 

3.60 

.96 


9 

3.26 

.87 


10 

3.93 

1 .om 


11 

'3.17 

1.69 
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Table 1.3. Orbit 331, 33*K to 43®N Interim Minns Nominal Mode 


Channel column 


6t6 V \ 

2 
3 

6#6 H I 

2 
3 

lOf 7 V I 

2 
3 

5 

6 
7 

10*7 H 1 

2 
3 

M 

5 

6 
7 

18.0 V I 

2 
3 
H 

5 

6 
7 

e 

9 

10 

11 


mean 

STD. 0£V. 

1.59 

.35 

i .87 

♦ '17 

1 .ss 

.29 . 

3.78 

,90 ■■ 

.20 

.29 

-.39 

.*11 

-.19 

.39 

-,9H 

.31 

1 .50 

. 5 1 

1.90 

.13 

2.02 

.98 

2.19 

.75 

2.39 

.60 

1 .SB 

-.51 

9,76 

.99 

3.31 

.52 

2.39 

.70 

2.01 

.56 

1 .97 

.62 

2.11 

.19 

2.21 

.92 

1 .96 

.91 

9.06 

1 .09 

2.80 

.89 

3.12 

1 .09 • 

2.86 

1.72 

3.06 

1.79 

2.91 

.95 

2.89 

.99 

3.18 

1 • 1 9 

2.86 

1 .06 

2.73 

.73 

3.67 

.89 
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Table 1.2. Orbit 331, 33"*N to 43®N Interim Minus Cross Mode (Continued) 


Channel column 


37«0 V I 

2 
3 
H 

5 

6 

7 

8 
9 

10 

11 

12 
13 

15 

16 

17 

18 

19 

20 
21 

' 22 

37#0 ri I 

2 
3 

- H 

5 

6 

7 

8 
9 

10 
1 1 
12 
13 
IH 

15 

16 

17 

18 

19 

20 
21 
22 


MEAK 

STO. OEV* 

2.94 

.63 

2.26 

.43 

2.52 

.48 

1.98 

.48 

1.91 

.44 

2.22 

.^1 

2.13 

♦ 42 

2.08 

.35 

2.Q8 

.41 

2*18 

.48 

2. 10 

.47 

1.97 

.43 

2.09 

.43 

1.93 

.43 

1.94 

• 38 

1.96 

.39 

1.90 

.38 

1.96 

.31 

1.81 

.41 

2.07 

.37 

2.08 

.40 

2.20 

.58 

3.5h 

.71 

3.7l 

.65 

3. 1 1 

.50 

3.69 

.57 

4.25 

.46. 

3.34 

.51 

3.49 

.54 

3.10 

.51 

3.17 

.52 

3.14 

.50 

2.90 

.58 

2.94 

.55 

2.93 

.53 

3.04 

.51 

2.99 

.49 

3.19 

.49 

2.87 

.54 

3.36 

.43 

3.61 

.69 

3.33 

1.43 

3.88 

1.12 

3.47 

1.90 
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Table 1.2* Orbit 331, 33®N to 43®N Interim Minus Cross Mode (Continued) 


Channel 

UtO H 
21 *0 V 
21 tO H 


column 


mean 


5TDf 


1 

2 

3 

5 

6 
7 
B 

9 

10 
1 1 

1 

2 

3 

5 

6 

7 

8 

9 

10 

I I 

1 

2 

3 

5 

6 
7 

e 

9 

10 

I I 


^^30 

.37 

3*98 

.2M 

3t IB 

• 28 

3# 10 

.27 

3.03 

.26 

3.00 

.23 

3.09 

.25 

3. 12 

.22 

3.21 

*25 

3.99 

.26 

3., 7 8 

.3h 

3.99 

,79 

3.58 

.58 

3.67 

,56 

3.61 

,H7 

3.98 

,H7 

3.S5 

.‘<8 

3.52 

,52 

3.62 


3.58 

,'<6 

3.59 

.52 

3f 33 

,5.7 

9.23 

.77 

3.50 

.60 

3.97 

.63 

3.70 

,5l 

3.37 

.55 

3.13 

.53 

3.27 

,5^ 

3.93 

.*18 

3.. 35 

, **8 

3.77 

.58 

3. 55 

.68 
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Table 1.2. Orbit 331» 33“N to 43°N Interim Minus Cross Mode 


channel 

column 

hean 

STD. dev* 

6*6 V 

1 

1 «S9 

.09 


2 

J » 45 

• 09 


3 

1.64 

• 0 6 



2.41 

. 17 

4t6 H 

1 

.20 

.06 


2 

-.00 

.08 


3 

.01 

.06 


H 

15 

.09 

50«7 V 

1 

1.50 

.07 


2 

1 .84 

.08 


3 

1.93 

.-08 


H 

2.05 

.08 


5 

2. 15 

.09 


6 

2.09 

• 08 


7 

2.73 

« 10 

10.7 H 

1 

3,12 

.07 


2 

2.39 

.07 


3 

2.09 

.06 



2.03 

.06 


5 

2.11 

.09 


6 

2.22 

.09 

• 

7 

2. 12 

.12 

I0»o V 

1 

2. 89 

.39 


2 

2.69 

.27 


3 

2.91 

.29 



2.78 

.25 


5 

2.84 

.25 


6 

2,73 

.24 


7 

2.82 

.27 


8 

2.84 

.22 


9 

2.79 

♦ 27 


10 

2.63 

.25 


1 1 

,2.79 

.34 
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Table l.l* Orbit 331, 33®N to 43*N Interim Minus Box Mode (Continued) 


Channel column 


37*0 V i 

2 
3 

5 

6 

7 

8 
9 

10 
I I 
12 
13 
\H 

15 

16 

17 

18 

19 

20 
21 
22 

37*0 H 1 

2 

3 

4 

5 

6 

7 

8 
9 

ID 
I 1 
12 
13 
IM 

15 

16 

17 

18 

19 

20 
21 
22 


mean 

STD. dev* 

2.73 

.90 

2.28 

.37 

2 »Hi 

• 3 6 

1 » 97 

.36 

1.83 

.35 

2.20 

• 36 

2.10 

.35 

2.06 

.29 

2. 16 

.31 

2.12 

• 39 

2.09 

.39 

1 .99 

.36 

2.03 

.33 

1.89 

.,30 

1 .93 

.28 

1.85 

.28 

2.00 

.28 

1 .9*» 

.27 

1.75 

.33 

2.05 

.30 

2.12 

.33 

2.07 

.33 

3.59 

•96 

3.67 

.95 

3.15 

.39 

3 . 75 

.92 

‘1.17 

.92 ■ 

3.37 

.92 

3.92 

.92 

3.20 

.36 

3.13 

.39 

3.11 

.92 

2. <J6 

. 90 

2.9H 

.90 

2.91 

.38 

3.05 

.37 

2.97 

.35 

3.16 

.36 

2.93 

.37 

3.32 

.36 

3.63 

.91 

3.92 

.52 

3.66 

.91 

3.73 

.78 
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Table 1.1. Orbit 331, 33“N to 43®N Interim Minus Box Mode (Continued) 


Channel 

ia»o H 


2l tO V 


21 .0 H 


COLUMN mean 5T0« DEV* 


1 

H .26 

.25 

2 

3 « 5 l 

.24 

3 

3 • 1 6 

• 24 


3 t 07 

*22 

5 

3*05 

.21 

6 

3 « 0 <( 

*20 

7 

3*00 

*21 

8 

3.11 

*20 

9 

3*22 

* 2 l 

10 

3*45 

*21 

1 1 

3*72 

• 18 

1 

3*70 

*58 

2 

3*60 

♦ 57 

3 

3*68 

* S 5 

H 

3.62 

.48 

5 

3*50 

■ .46 

6 

3*58 

.47 

7 

3 • 55 

*48 

8 

3 * 60 

• .4m 

9 

3 * 60 

.46 

iO 

3.50 

*49 

1 1 


.47 

I 

*».23 

* 66 

2 

3*55 

. 6 1 

3 

3.93 

.60 


3.68 

• .53 

5 

3.38 

* 54 

& 

3 * 15 

. S -2 

7 

3.26 

.52 

8 

3.41 

.48 

9 

3.36 

.48 

10 

3.76 

♦ S 3 

1 1 

3*67 

*46 
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Table Orbit 331, 33®N to 43^N Interim Minus Box Mode 


Channel coluhn 


V . 1 

2 
3 

6 H I 

2 
3 
H 

1 0 # 7 V I 

2 
3 

5 

6 
7 

I0t7 H 1 

2 
3 

5 

6 
7 

18»0 V 1 

2 
3 
H 

5 

6 

7 

8 
9 

10 

11 


mean stq# dev* 


1*60 

.06 

I *22 

♦ 07 

1 »30 

.06 

1 • S 7 

-.12 

.21 

.05 

.19 

.06 

.23 

.09 

• 3 6 

.05 

1 .52 

.09 

1 .80 

.05 

1.90 

.09 

1 .99 

.0,9 

2.08 

.05 

2.02 

.05 

2.-06 

.07 

3.03 

.09 

2*HH 

.09 

2.19 

♦ 09 

2.08 

.09 

2.09 

.05 

2.25 

.05 

2.35 

.06 

2.51 

.26 

2 .70 

.26 

2.93 

. 2 *6 

2.83 

.23 

2.89 

.23 

2.75 

.22 

2.85 

.23 

2.89 

.22 

2.83 

.29 

2.69 

.29 

2*66 

*25 
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truth data sets. These data sets can then be used along with 
geophysical models and appropriate assumptions about atmospheric water 
content to yield a better determination of the bias for each SMMR 
channel. To facilitate this process, it is recommended that, as new 
APC output is produced, the LOCATE program be rim in a production mode 
to accumulate additional buoy hits. 

6.0 NEW TECHNOLOGY 

No new technology has been developed in the course of this study. 

7.0 REFERENCES 

1. Kitzis, S.N. and Kitzis, J.L. , ’’Evaluation and Analysis of 
SEASAT-A SMMR APC Algorithm; 6 . 6 - GHz T^ vs. Te„-rforp 
Comparison Results”, March 16, 1979. 

2. Kitzis, S.N. and Kitzis, J.L., ’’Evaluation and Analysis of 
SEASAT-A SMMR APC Algorithm: T^ Measured vs. Tg Calculated 
Comparison Results”, April 13, 1979. 

3. Kitzis, J.L. and Kitzis, S.N., ’’Evaluation and Analysis of 
SEASAT-A SMMR APC Algorithm: Interim T 3 vs. Final with 
Gaussian Coefficients”, May 25, 1979. 

4. Kitzis, J.L. and* Kitzis; S.N., "Evaluation andr Analysis of 
SEASAT-A SMMR APC Algorithm: Interim Mode T^ vs. Cross and 
Nominal Mode Tg”, July 27, 1979. 


15 



(2) Initial analysis carried out by H. Hofer of data produced by the 
software described in section 3.2 indicates that there are still 
significant inter-channel biases in the SMMR data (i.e. residual 
biases between SMMR Tg values and model-predicted brightness 
temperatures from surface truth) . Simulation runs indicate that 
these inter-channel biases are not introduced by the APC 
algorithm. Therefore, these biases are either the result of 
errors in the model-predicted brightness temperatures or else 
biases in the data which is input to the APC. 

(3) Over open ocean, the only differences observed between the four 
APC modes are offsets which vary with the mean level. At 
low T 3 levels, the interim mode values are above the other three 
modes, while at high Tg levels the reverse is true. This chang- 
ing 'offset varies smoothly with increasing brightness 
temperatures ♦ 

(4) The box mode outputs are generally similar to those of the interim 
mode, except that the box mode removes some of the sidelobe con- 
tributions from nearby land. No "ringing" effects are observed 

in the box mode values. Thus, the box mode is superior to the 
interim mode near land without suffering any ill effects from 
"ringing." 

(5) The cross mode appears to remove more of the sidelobe contributions 
near land than does the box mode. TJnfortimately, the cross mode 
exhibits "ringing" effects for some channels, thereby making it 
less desirable than the box mode. 

( 6 ) The nominal mode still shows severe "ringing" effects for all 
channels. Although this mode shows the potential for signi- 
ficantly better resolution than the other modes, the "ringing" 
renders it unreliable at the present time, 

5.0 RECOMMENDATIONS 

All of the previous questions regarding the SMMR brightness tempera- 
tures have been answered with the possible exception of two. One 
involves the "ringing" effects observed in the nominal and, _to_a much 
lesser degree, in the cross modes near land. _At_ present, this ^ 

"ringing" phenomenon is not completely understood and should be further 
investigated unless a final decision is made to completely disregard 
the nominal and cross modes. This effect can best be studied by simu- 
lating a brightness temperature scene which is half ocean and half land, 
which requires a more sophisticated simulation approach than that 
currently in use. 

The second remaining question concerns the inter-channel Tg biases. 

These biases should be studied by using the software described in sec- 
tion 3.2 to assemble a large number of matched spacecraft and surface 
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@ASG,A JK*ENI. 

@ASG,A JK*TABLE. 

@ASG,A 1135*WINDFIErD. 
euSE 10.,1135*WINDFIEID. 

(?ASG,T TAPE. ,U9H,Sxsx3cR (Tape containing Tj data) 

(3REWIND TAPE. 

@USE SMMRTAPE., TAPE. 

(§CAT,P ODTPUTl. (Output file for T^ data) 

®ASG,A OUTPUTl. 

@USE 23., OUTPUTl. 

@CAT,P OUTPUT2. (Output file for surface truth data) 

@ASG,A OUTPUT2. 

@USE ll.,0UTPUT2. 

@XQT JK*ENI. CORAL 

stimes - 

LWIND = .TRUE., LVAPOR = .TRUE., 

$END 

@ADD JK*TABLE. 

The output file of Tg data contains two card images corresponding to 
each grid-1 cell location. The first card image contains an ID number 
for the point and the ten SMMR brightness temperature measurements. 

The format of this card is 15, 10F7.2 . 

The second card image- contains the same ID nirmber and the ten grid 
column numbers corresponding to the Tg measurements. This card has the 
format 15, 1017. The ID numbers range from 101 through 116 for the 16 
points from the first block of SMMR data processed, then from 201 
through 216 for the points from the second block of data, etc. 

The output file of surface truth data contains a single- card image 
corresponding to each grid-1 cell location. This card contains the 
ID number for the point, the sea surface temperature in degrees 
Celsius, the water vapor content in grams /cm^, and -the wind speed in 
meters /second. The format used is 15, 3F7.2 . 

4.0 CONCLUSIONS 

Several conclusions may be drawn from the results discussed in 
section 3.0: 

(1) ’Most of the problems in the APC outputs discussed in previous 
reports have now been eliminated or reduced to insignificant 
levels. These problems all involved unrealistic variations in 
brightness temperatures across the SMMR swath. These include 
uhmodeled cos 0 effects, cross-track gradients induced by the 
sidelobe correction procedure, and large biases at the edges of 
the swath also caused by sidelobe corrections. These problems 
are no longer apparent in the brightness temperature output of 
the APC. 
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choose to have the temperatures for all channels taken from grid 1 , or 
alternatively to have temperatures for each channel taken from those 
cells of its best grid which lie closest to the grid*-l cell locations. 

In either case, the grid 1 locations are used to interpolate the surface 
truth fields to produce corresponding surface truth measurements. If 
a given location is outside the limits of a field, the corresponding 
surface truth parameter is set to zero. 

The namelist $TIMES is required to execute the program and contains the 
following variables: 

Variable 


. Name 

Type 

Dim. 

Default ' 

Description 

ITl 

' Int. 

4 

0 . 0 , 0,0 

Beginning time (Days ,Hrs ,Min , Sec) 

IT2 

Int. 

4 

366,0,0,0 

Ending time (Days,Hrs, Min, Sec) 

ALLGRl 

Log. 

1 

.FALSE. 

♦TRUE, if T 3 measurements for all 
channels are to be taken from grid 1 . 

LWIND 

Log. 

1 

.FALSE. 

♦TRUE, if wind speeds are to be inter- 
polated from a wind field. 

LVAPOR 

Log. 

1 

.FALSE. 

.TRUE, if water vapor values are to be 


calculated . 


The namelist $TABLE contains the sea surface temperature field data and 
is made available to the program by 0ADD^ing the file JK^TABLE to the 
runstream. In addition to printed output, the program produces two 
output files, one containing T^ data and one containing surface truth 
data . 

The following is a sample runstream for orbit 1135, assuming that Tg 
data is to be read directly from a nine-track tape. 

0ASG,T TAPE.,T,Q160K 
0EZIO,RS JK*ENI. 

0REWIND TAPE. 

0FREE TAPE. 

0ASG,T TAPE. ,T,Y317R 
0EZIO,RS JK^TABLE. 

0REWIND TAPE. 

0FREE TAPE. 

0ASG,T TAPE.,T,W498R 

0EZIO,RS TAPE. , 1135*WINDFIELD . , 1298*WINDFIELD . 

0REWIND TAPE. 

0FREE TAPE. 
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these are not measured by buoys and are simply set to zero* The fields 
are as follows: 


Field # 

Parameter 

Units 

Format 

1 

Air temperature 

Tenths of degrees Celsius 

. 14 

2 

Dew point 

Tenths of degrees Celsius 

14 

3 

Barometric pressure 

Tenths of millibars 

15 

4 

Wind speed 

Hundredths of meters /sec. 

14 

5 

Wind direction 

Tenths of degrees from* north 

14 

6 

Weather code 

— 

11 

7 

Visibility 

Tenths of nautical miles 

13 

8 

Precipitation 

Millimeters 

14 

9 

Solar radiation 

Hundredths of Langleys /min 
(wave lengths from .6 to 
4.0 microns) 

13 

10 

Solar radiation 

Hundredths of Langleys /min 
(wave lengths from 4.0 to 
50 microns) 

13 

11 

Significant wave height 

Tenths of meters 

13. 

12 

Average wave period 

Tenths of seconds 

13“ 

13 

Average wave direction 

Degrees from north 

13 

14 

Highest crest 

Tenths of meters 

13 

15 

Deepest trough 

Tenths of meters 

13 

16 

Sea ’ surface temp * 

Hundredths of degrees Celsius 

14 

17 

Salinity 

Thousandths of parts per 
thousand 

15 

18 

Conductivity 

Thousands of millimhos/cm 

15 


The third and fourth cards respectively contain the maximum and minimum 
values measured during the six-hour window* The fifth card contains the 
standard deviations of the measurements about their mean values* These 
three card images use the same format as the second card* 


3*2.3 CORAL User^s Guide * The program CORAL uses sea surface temperature 

fields, wind fields, and a simple water vapor model to pt'o^uce’ corres- 
ponding sets of surface truth values and SMMR brightness temperature 
measurements. The current sea surface temperature field covers a 
region in the northeastern Pacific between 21 degrees and 55 degrees 
north latitude, and between 179 degrees and 245 degrees east longitude. 
It represents the NMFS monthly average temperatures for September 1978. 
Figure 15 shows the full extent of the sea surface temperature field. 
Wind field values have been produced by Vince Cardone and are avail- 
able for SEAS AT orbits 1135 and 1298. Values for water vapor content 
are calculated as a fiinction of latitude. 

Sixteen brightness temperature measurements for each channel are 
extracted from each block of SMMR data,. These* measurements are asso- 
ciated with the centers of the sixteen grid-1 cells* The user may 
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In addition to printed output, a FORTRAN formatted file is produced • 

An example of the printed output appears in figure 14* The nonzero 
measurements are respectively air temperature (C) , pressure (mb) , 
wind speed (meters/sec), wind direction (degrees from north), and sea 
sxnrface temperature (C)* The value of 1 for JSTAT indicates that 
measurements are available both before and after the time of overflight 
within the 6-hour window centered on overflight • A value of 2 would 
indicate that measurements were found on only one side of the over- 
flight time. The value of KTOT indicates the number of sets of mea- 
surements found within the 6-hour window centered on overflight. When 
no measurements are available, an appropriate message is printed. 

The following is a sample runstream: 

0ASG,T TAPE.,T,Q160R 

0EZIO , RS TAPE . , JK*ENI . , SEASAT-ADF*BFILE . 

0REWIND TAPE. 

0FREE TAPE. 

0USE BFILE.,SEASAT-ADF*BFILE. 

@ASG,A JK^ENI. 

@ASG,A BFILE. (This buoy file must be assigned to the run under the 

name BFILE.) _ 

(5ASG,A pmPUT.^_ (Output fde)_/ 

0USE 22". , OTJTPUtV ' ^ ' 

0ASG,A HITFILE. (File produced by LOCATE) 

@USE 23., HITFILE. 

0XQT JR*ENX. HATCH 

The output file contains either one or five card images corresponding 
to each overflight time encountered in the hit file from LOCATE, 
depending on whether or not buoy measurements are available within 
3 hours of overflight. The first card image is formatted as follows: 

‘(Buoy ID) (Days,Hrs,Min,Sec) JSTAT NTOT 

110 5X 415 5X 5X II 

A value of 3 for JSTAT indicates that buoy measurements are not avail- 
able and that no additional cards have been produced for this over- 
flight time. A value of 1 or 2 indicates that measurements are ^ 
available, and in this case, the 4 cards immediately following contain 
information regarding this hit. 

The second card image contains the interpolated measurements. There 
are fields for eighteen separate physical parameters, but many of 
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The following is a sample mnstream, when reading Tg data from a 
9-track tape* 

@ASG,T TAPE.,T,qi60R 
@E2I0,RS TAPE*,JK*EN1* 

@RERIND TAPE. 

@FREE TAPE. 

@ASG,A JR^ENI. 

@ASG,A OLBHITS. (file of previously accumulated hits) 

@USE 25.>OLDEITS. 

@CAT,P HITPILE. (output file - will include old and new hits) 

@ASG,A HITFILE. 

@USE 23., HITFILE. 

@ASG,T TAPE.,U9H, SxxxxR (Tape containing Tg data) 

@REWINI) TAPE. 

@USE SMMRTAPE., TAPE. 

@XQT JK^ENI. LOCATE 
$TIMES 

MERGE — • TRUE » , 

$ENI) 

Note: If MERGE - .FALSE., no file of previously stored hits is -- 

required. 

The output file contains 3 card images corresponding to each hit 
encountered. The first card is formatted -as follows: 

(Buoy ID) (Days, Hrs., Min., Sec.) (Seconds from beginning of year^ 
11^ 5X 415 ' 5X FIO.O 


The second card image contains the ten brightness temperature mea- 
surements and uses the format 10F7.2. The third card image contains 
the ten grid column numbers, and uses the format 1017. 

3.2.2 _MATCH User^s Guide . The program MATCH outputs buoy surface truth data 
corresponding to overflight times listed in a **hit” file produced by 
LOCATE. Surface truth values are extracted from a buoy file covering 
the period between July 1 and October 31, 1978, and are interpolated to 
the exact time of overflight whenever possible. Measurements must be 
available within 3 hours of overflight in order to produce surface 
truth ‘output. If measurements^ are available both before and after the 
time of overfligh t, a linear interpolat ion is performed. If measure- 
ments are only available on one side of the time of overflight, those 
closest to the time are used. The program also outputs the maximum 
and minimimi values measured within 3 hours of overflight and the 
standard deviations of the measurements about their mean values during 
the same period. Measurements contained in the buoy file include air 
temperature, air pressure, wind speed and direction, sea surface tem- 
perature, and sometimes significant wave height and average wave period 
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to yield valid statistical results. In order to automate the process 
of accumulating large sets of spacecraft and surface truth data which 
are co-located in time and space, we have developed a package of three 
programs. These programs are: 

(1) LOCATE - This program identifies times at which the SMMR swath 
passes over the locations of a predetermined set of buoys. 

(2) MATCH - This program produces buoy measurements corresponding 
to the overflight times output by LOCATE. 

(3) CORAL “ This program uses digitized surface truth fields 

to produce corresponding sets of surface truth measurements and 
SMMR data. 

These programs have been designed with the goal of minimizing required 
user preparation and involvement. They may be used in a production 
environment in which the only inputs which change for each run are the 
file names/tape II)*s to be processed. Since it is expected that a 
large number of users may be interested in this software, a user’s 
guide for each program is provided in the following sections. This 
software has been thoroughly tested and has been used to process all 
the SMMR T 3 data produced for the SMMR Mini-II workshop. These pro- 
grams are easy to set up and are fast-running. 

3.2.1 LOCATE User^s Guide . The program LOCATE identifies ’times at which the 
SMMR swath passes over the location of any of .the buoys listed in 
Table 3. The outputs of the program include buoy ID, time of over- 
flight, SMMR brightness temperature measurements for all 10 channels, 
and the grid column numbers corresponding to these measurements for each 
buoy "hit" encountered. The Tg measurements for each channel are taken 
from that channel’s best grid and represent the grid cell closest to 
the location of the buoy. A buoy falling outside the SMMR swath may 
still result in a hit if the buoy is within 13.5 kilometers (one half 
of a grid 4 cell-width) of the swath’s edge. In addition to printed 
output, the program produces a FORTRAN formatted file for use in a 
subsequent program. 


A namelist input. STIMES is required to execute the program, and includes 
the following variables: 


Variable 

Name 

Tjrpe 

Dim. 

Default 

Description 

ITl 

Int. 

4 

0 , 0 , 0,0 

Starting time for processing Tg data 
(days, hours, minutes, seconds) 

IT2 

Int. 

4 

365,23,59,59 

Ending time for processing Tg data 
(days, hours, minutes , seconds) 

MERGE 

Log. 

1 

.TRUE. 

.TRUE.’ to merge all hits encountered 


with file of previously accumulated 
hits. 
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Observations regarding these plots are listed below: 

(1) The effects of the adjacent land on the interim mode data can be 
seen as a rise in T 3 values occurring between 38® and 40® north 
latitude. The horizontal channels show larger effects than the 
vertical channels. In fact, the interim 18V and 21V channels 
(figures 11.5 and 11.7) show hardly any response to the nearby 
land. This behavior is unchanged from that previously observed 
in references 3 and 4. 

(2) For all channels except 6.6 GHz, the only difference between the 
box and interim modes is a constant offset. For the 6.6 GHz 
channels, this offset diminishes slightly as the SMtIR swath 
approaches the point of closest proximity to land. These obser*- 
vations imply that the box mode sidelobe corrections are not 
improving upon the performance of the interim mode in this par- 
ticular case. In fact, the observation concerning the 6.6 GHz 
channels seems to suggest that the box mode may be more sensitive 
to the nearby land than the interim mode. 

(3) The behavior of the cross mode is somewhat similar to that of the 
box in that the cross mode values differ from the interim mode., 
values by a nearly constant offset. However, this offset 
increases near the point of closest proximity to land for the 
6.6 H and 18 H channels. This implies that the cross mode is 
less sensitive to the nearby land than the box or interim mode 
for these channels. Unfortunately, the cross mode data for the 

10 . 7 H channel exhibits some over- and under- compensate on which, is 
reminiscent of the ’^ringing*’ discussed for the coastline crossing 
case, 

. (4) As noted in prior reports, the nominal mode data exhibits more*‘ 
noise than the other modes. In addition, the nominal mode shows 
severe over- and under-compensation effects near land for all the 
horizontal channels. This is particularly evident for the 6.6 H 
and 10.7 H channels (figures 13.2 and 13.4). .This severe over- 
and under-compensation is undoubtedly related to the "ringing” 
effects previously discussed. 

3.2 Creation of Matched Surface Truth and Spacecraft Data Sets . Refer- 

ences 2 and 3 document attempts to determine instrument biases in the 
SMMR data by comparing SMMR Tg values with surface truth measurements. 
Such comparisons are difficult to perform due to the frequent lack of 
adequate atmospheric water measurements. This problem may be partially 
alleviated by first assembling large quantities of surface truth data. 
Measurements corresponding to clear atmospheric conditions can then 
be extracted from the large ensemble, or alternatively, mean values of 
atmospheric water content can be assumed to hold for the complete data 
set. In either case, a sufficiently large data set must be available 
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Observations regarding the complete set of data are listed below: 

(1) The box, cross, and nominal mode values are generally below the 
interim values over the ocean. This is not the case over land, 
where the relationship is reversed for the lower frequency 
channels . 

(2) As expected, the transition from ocean Ho land is sharpest for the 
37 GHz channels and becomes less sharp xTxth decreasing frequency. 
In addition, the ocean data exhibits more noise for the higher 
frequency channels than for the lower frequency channels. This 

is due to the better resolution of small local weather features 
afforded by the higher frequencies. 

(3) ’^Ringing" is clearly observed in^ the nominal m data near the 
coastline crossing (figures 10.1 through 10.10). It is present 
to a much lesser extent in the cross mode data (for 10. 7H and 
18V), and is not present at all in the box and interim modes. 

(4) The differences observed between the interim mode and the other 
three modes generally increase as the distance from land dimin- 
ishes. This increase in the differences is due to a gradual rise 
in the interim mode values as the swath approaches land while the 
box and cross modes maintain more constant levels. This is good 
evidence that the side lobe correction procedure is truly removing 
contamination from nearby land. 

(5) Table 2 summarizes the differences observed for orbit 1212, over 
open ocean and adjacent to the coast. Note that the land con- 
tamination in the interim data extends farthest from the coast 
for the lower frequency channels (e.g. 7 grid-1 cells 1000 km. 
tqv 6.6 GHz versus 4 grid-4 cells ^ 100 km for 37 GHz), which is 
3 ^^^l^cti^ o^ their wider beam widths. Also, note that the 
cross mode appe ars t o removF mor e~^f ' the”Tahd c on t agination 
than does the box mode, whichT is evidenced by the large’r’ dif-’ 
ferences adjacent to land for the cross mode. Unfortunately, 
the cross mode appears to over-correct for land (i.e. "ring”) 

in the 10. 7H and 18V channels. This over-correction becomes 
much more severe in the nominal mode for all channels. 

3.1.3 Parailellng Coastline Case . We have analyzed data from orbit 1212 
near 38° north latitude as an example of the SMMR swath following a 
course parallel and extremely close to land. The right-hand edge of 
the swath runs parallel to the California coastline and comes within 
50 kilometers of San Francisco. There are two points of closest 
approach to land, one at 38.5° and the other at 40° north latitude. 
Figures 11, 12, and 13 contain plots of Tb values versus latitude for 
the four APC modes . 
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(2) Table 1 also shows that the interim mode values are consistently 

higher than the corresponding box, cross, and nominal mode values 
for all channels • The mean differences tend to increase with 
increasing frequency, and vary between 0® and These differ- 

ences are approximately the same for the box, cross, and nominal 
modes » 

(3) The outer cells for the 18, 21, and 37 GHz channels tend to be 
somewhat lower than the center cells in the box, cross, and 
nominal mode outputs. 

(4) The cross-track gradients observed in the T 3 data appear to be 
very similar for all four APC modes. A slight enhancement of the 
interim mode gradients can be seen in the box, cross, and nominal 
modes (see figures 4 through 7) . These gradients for all modes 
appear to be geophysically realistic, and thus indicate that the 
final cos 3 estimates have been properly implemented. 

(5) 'The box, cross, and nominal modes appear to introduce small 

(<1®K) opposing cross-track variations in the 10.7 V and H data 
(see Table 1) . This is the only remaining evidence that the 
sidelobe correction procedure is introducing undesired effects. 
Note, however, that ail of the previously observed problems 
(see reference 4) have been corrected. 

( 6 ) As noted previously.. in references 3 and 4, the -difference between 
the interim mode and the other modes decreases with an increasing 
average Tg level. In fact, if the Tg level increases high enough, 
the box, cross, and nominal mode values will be consistently 
above the interim mode data. An excellent example can be seen in 
figures 1 , 2 , and- 3, where data is plotted through the eye of 

hurric an e Fico. N ote t hat t he in terim values are abo ve t he other . 

mode values away from Fico, but are. below them near the eye of 
Fico, with a smooth transition in between. 


(7) An additional comment is in order regarding the pass over Fico. 
The nominal mode appears to achieve a better resolution of the 
eye of the hurricane than do the other modes. ‘This is readily 
apparent in figures 3.1 through 3.10, and especially in fig- 
ures 3.3, 3.4, and 3.6. The box and cross modes generally show 
the eye to the same extent as the interim mode, whereas the nomi- 
nal mode shows a significant improvement in resolution. 

3.1.2 Coastline Crossing Case . *We have analyzed data from orbits 1178 and 
1212 as examples of the SMMR swath crossing a coastline. We have 
produced plots of Tb versus latitude corresponding to these two orbits 
for different positions across the scan and for each of the APC modes*. 

A set of expanded-scale plots produced for the left-hand edge of the 
SMMR swath for orbit 1212 are displayed as figures 8 , 9, and 10. The 
sharp rise in Tb values near 50 is due to the SMMR swath crossing 
the Kenai Peninsula on the Alaskan coast. 
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A third analysis tool calculates mean Tg values for each cell of the ten 
channels for two modes of the APC, and associated standard deviations^ 

The program, also calculates the mean differences between the two modes 
and associated standard deviations. This program is typically used to 
compare interim mode results with results of one of the other modes. An 
example of the program’s output may be found in Table 1. 

The comparisons presented here have been performed for the purpose of 
evaluating the effects of the different sidelobe correction procedures. 

The following results are expected from these types of comparisons: 

(1) All four modes of- the APC should produce generaily similar results 
over open ocean, in terms of both mean Tg levels and any gradients 
which may exist in the data. However, the nominal mode is expected 
to exhibit more point-to-point variation than the interim mode, 
with the box and cross modes being intermediate between the two. 

(2) When crossing a land-sea interface, the sharpest transition is 
expected for the nominal mode, a somewhat less sharp transition for 
the cross and box modes and the smoothest transition for the interim 
mode. The cross and box modes are expected to remove the ringing 
phenomena observed in the nominal mode and discussed in references 

3 and 4. 

(3) When paralleling a coastline, the nominal mode should be least sen- 
sitive, the cross and box modes slightly more sensitive, -and the 
interim mode most sensitive to the nearby land. 

This section contains the results of a direct comparison of interim mode 
data with nominal, cross, and box mode data from the APC. Separate com- 
parisons are made for the open ocean, land-crossing, and land-paralleling 
•cases listed above. 

3.1.1 Open Ocean Case . Por the open ocean case, we have analyzed data from two 
SEASAT passes over the Gulf of Alaska, orbits 1178 and 1212, and a third 
pass over hurricane Pico near Hawaii, orbit 331. Por each pass, we have 
produced plots of brightness temperature versus latitude and of cross- 
track gradients versus latitude. In addition, we have calculated mean 
differences between the various modes for each grid^ cell. Examples of 
these three types of output for the hurricane pass may be found in Pig- 
ures 1 through 7 and Table 1. 

Observations regarding the complete set of data are listed below: 


(1) As expected, the nominal mode data exhibits more noise than any of 

the other modes. As seen in Tables 1.1, 1.2, and 1.3, the noisiness 
of the box and cross modes is only slightly greater than that of the 
interim, while the nominal exhibits significantly more no5vSe than 
the others. 
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Report 4 - The fourth report (reference 4) extended the evaluation to the 
Nominal and Cross modes of the APC. Both of these modes were compared 
directly with Interim APC output for the three environmental regimes. As 
observed for the Gaussian APC, the Nominal and Cross mode Tb values were 
biased low with respect to Interim values a and exhibited cross-track 
gradients significantly different from those observed in Interim 
output. In general, these effects were less severe for the Cross than 
for the Nominal APC. "Ringing^* was again observed in Nominal output near 
coastline crossings. Two parallel efforts showed that the world T 3 map 
agrees well with actual SMMR data, and that the cos 3 correction effec- 
tively removes the previously observed opposing cross-track gradients. 

Report 5 - This present report evaluates the latest versions of the 
Interim, Cross, Box, and Nominal modes of the APC. Most of the problems 
observed in report 4 have been corrected. The biases between the Interim 
and the other modes are significantly reduced, and the gradients pre- 
viously introduced by the sidelobe correction procedure are no longer 
apparent. ’’Ringing" is still observed in the Nominal, is greatly 
reduced in the Cross, and does not appear at all in the Box APC. For 
the first time, cases have been observed which show that the implemented 
sidelobe correction procedure removes land contamination effects such 
that the Nominal, Cross and .Box mode data is preferable to the Interim 
data when near land. In addition, due to the need for larger surface 

truth data sets, new software has_been_developed which facilitates__the 

* production of matched space_cr^t m'd surface truth data _ sets _ which are 
co-located in time and spacer. ...f. 

3.0 TECHNICAL DISCUSSION 

3.1 Inter comparison of APC Modes . The comparisons presented here are a - 
continuation of the work previously discussed in section 3,1 of refer- ^ 
ence 4. Several software analysis tools have been used to facilitate 
these comparisons. 

One of these tools produces plots of brightness temperatures versus 
latitude for two modes of the APC. Typically, these plots display interim 
mode data along with either nominal, cross, or box mode data. Examples 
• of these plots are shown in figures 1.1 through 1.10. For all plots pre- 
sented here, the s 3 nnbol "I" represents interim mode data, while "F" 
represents either box, cross, or nominal mode data, as indicated in the 
plot title. 

A second analysis tool produces plots of T^ cross-track gradients versus 
latitude. These gradients are calculated by fitting a straight line to 
each row of Tg data cells. The gradients are plotted in units of degrees 
Kelvin per grid cell. Examples of these plots are shown in figures 4.,1 
through 4.5. In each plot, the symbol "V" represents gradients for the 
vertical polarization and "H" represents gradients for the horizontal 
polarization for the frequency indicated in the plot title. 
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Of course, each of these versions has incorporated changes over time as 
problems have been uncovered and the resulting corrections implemented ♦ 

Por example, the Nominal, Cross, and Interim modes of the APC discussed 
in this report are slightly different from those of the same names dis- 
cussed in the last report (reference 4) . These changes are a reflection 
of problems identified in this prior report, as well as problems 
mcovered in independent investigations* Even though the present report 
concludes this particular evaluation effort, it is expected that*' the APC 
will continue to change following the publication of the report There- 
fore, the results presented here should not be construed as representing 
the final state of the APC algorithm, but only its status as of this time. 

The following is a short history of the valiant attempts of this evalua- 
tion to keep pace with the ever-changing APC algorithms* 

Report 1 - This first report (reference 1) documented the initial evalua- 
tion of the Interim APC 6*6 GHz channels* SMMR data was compared with 
sea surface temperature values obtained from the NMFS monthly average sea 
surface temperature map" for September, 1978. This comparison yielded the 
first indication of the existence of opposing cross-track gradients. A 
brief survey of data suggested that the source of the gradients may 
lie in the T^ data. 

Report 2 - The second report (reference 2) extended the evaluation to all 
ten of the Interim APC channels. This portion of the evaluation estab- 
lished the basic software analysis tools to be used throughout the fol- 
•lowing stages of the study. SMMR data was compared with NBBO buoy 
measurements, using a geophysical model to convert the buoy measurements 
into predicted brightness temperatures. This initial attempt at deter- 
mining Tg biases was frustrated by a lack of atmospheric water vapor and 
cloud liquid water measurements. However, this comparison still indicated 
the existence of Tg biases which vary from channel to channel, with the 
horizontal channel biases being more positive than the corresponding 
vertical channel biases. In addition, opposing cross-track gradients 
were observed for the 18 and 21 GHz channels as well as the 6*6 GHz 
channels . * . - 

Report 3 - The third report (reference 3) represented the first evalua- 
tion of the sidelobe corrections implemented in the Gaussian APC. 

Gaussian APC output was directly compared with corresponding Interim APC 
output for the three basic environmental regimes of open ocean, crossing 
a coastline, and paralleling a coastline. The Gaussian T^ values were 
found to be significantly lower than corresponding Interim values, and 
in addition, gradients observed in the Gaussian data were significantly ^ 
different from those observed in Interim data. Both observations sug- 
gested that the sidelobe correction procedure was introducing errors 
into the data. This conclusion was reinforced by the first observa- 
_tioM of ^'ringing, ^en in the Gaussian T 3 data at coastline crossings. 

In a parallel effort ~first”estimateV'bf "cos gWalue^ were "obtained "by^ 
fitting T^ data to the narrow beam antenna equations. These cos g 
values correlated well with the opposing cross-track gradients pre- 
viously observed in the Interim APC. 
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1.0 


SUMMARY 


This report contains an evaluation of the latest versions of the SEASAT-A 
SMMR Antenna Pattern Correction (APC) algorithm, and concludes the APC 
evaluation previously documented in references 1 through 4. The report 
focuses on two principal efforts: 1) the intercomparison of the interim, 
hox, cross, and nominal APC modes, and 2) the development of new software 
to facilitate the creation of matched spacecraft and surface truth data 
sets which are co-located in time and space. 

The more important conclusions of the report are: 

(1) Most of the problems in the APC outputs discussed in the previous 
reports have now been eliminated or reduced to insignificant levels. 

(2) The box mode is superior to the interim mode in that it removes 
some of the sidelobe contamination effects due to nearby land with- 
out exhibiting *Vinging." 

(3) The cross and nominal modes are unreliable near land due to the 
presence of "ringing” effects. 

2.0 INTRODUCTION 

The work described in this report concludes the current evaluation ’of the 
APC algorithm which is documented in references 1 through 4. Over the 
.course of this evaluation, the APC has evolved through several stages of 
development. The following presents a synopsis of this stepwise 
■evolution. 

a. The first version of the APC to be implemented is known as the ^ 
Interim APC. This version of the algorithm does not include correc- 
tions for Earth sidelobe contributions, although it does correct for 
sidelobe contributions from cold space-. 

b. The Gaussian APC is the first version which applies Earth sidelobe 
corrections. The actual antenna pattezms are not used, but are 
approximated by smooth Gaussian functions. 

c. The Nominal APC is identical to the Gaussian version except that the 
actual antenna patterns are used in performing Earth sidelobe correc- 
tions. Corrections are applied for all cells outside the cell being 
corrected. 

d. The Cross APC is a refinement of the Nominal, developed for the 
purpose of eliminating algorithm sensitivity to small errors in the 
sidelobe contributions from the four cells nearest the one being 
corrected. 

e. The Box APC is a further refinement of the Nominal, which el iTm' -* 
nates sensitivity to small errors in all eight adjacent cells. 
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ABSTBACT 


This report contains an evaluation of the latest versions of the SEASAT-A SMMR 
Antenna Pattern Correction (APC) algorithm, and concludes the current APC evalua*- 
tion efforts The report focuses on two principal efforts: 1) the intercomparison 

of the interim, box, cross, and nominal APC modes, and 2) the development of new 
software to facilitate the creation of matched spacecraft and surface truth data 
sets which are co-located in time and space. The major conclusion of the report 
is that most of the problems discovered in earlier versions of the APC have now 
been corrected. 
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